1 



m 
o 
o 

(N 
CD 
0^ 



Nonperturbative Renormalization for Domain Wall Fermions and the 
Chiral Condensate 

Azusa Yamaguchi ^ 

'^Department of Physics, Columbia University, New York, USA 

We study the chiral condensate, {tpip), and various quark bilinear vertex functions for domain wall fermions at 
different lattice scales, with both the Wilson and DBW2 gauge actions, in both quenched and dynamical fermion 
simulations. We use the vertex functions to calculate renormalization factors within a non-perturbative scheme. 
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1. Introduction 

The RBC collaboration has performed simu- 
lations with domain wall fermions (DWF) and 
the DBW2 gauge action and presented results 
for weak matrix elements. This report focuses 
on the chiral condensate (ipip) and exploits non- 
perturbative renormalization calculations to get 
physical values for both quenched and dynamical 
fermions. The simulation details are given the Ta- 
ble Two lattice volumes are used, 24'' x 48 for 
= 3GeV, and 16'^ x 32 for the others. Domain 
wall fermion are used for all simulations with 
three values of Lg'- 16 for the Wilson gauge action 
and for the DBW2 gauge action with = 2GeV 
and l.BGcV, 10 for = 3GeV and 12 for dy- 
namical fermions. Among these, data and sim- 
ulation details for the Wilson gauge action and 
DBW2 gauge action with l.SGeV and 2GeV have 
already reported|l|. The DBW2, 3GeV quenched 
and dynamical fermion simulations with a^^ ~ 
l.SGeV are reported in this meeting by RBC|2j- 

2. {'ipip) measurement 

Using DWF, we can compute the value of 
the chiral condensate directly by evaluating 
(ipip) a series of quark mass and extrap- 
—rrires (see the section 4). 

= SGeV case, we use 
0.008, 0.016, 0.024, 0.032, 0.04 
10~^|31. For the dynamical 



olating to mf = 
For the quenched, 
5 fermion masses: 
and rures — 9.72 

case, three dynamical fermion masses 
0.02,0.03 and 0.04, and 5 valence masse 
0.01, 0.02, 0.03, 0.04 and 0.05 are used. In the dy- 
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Table 1 

Simulation details. Here "Q" means quenched 
while "D" means dynamical. The 3GeV simu- 
lation uses a 24-^ x 48 volume while a 16"^ x 32 
volume is used for the others. 
Gauge action Fermions Ls 



^(GeV) 



Wilson 
DBW2 
DBW2 
DBW2 
DBW2 



namical case, in order to get a physical value, 
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16 


1.922(40) 
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16 


1.97(4) 
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16 


1.31(4) 
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10 


2.86(9) 


D 


12 


1.80(7) 



val 



dyn 



we extrapolate along the line 
— TOics- In the dynamical simulation, rrires = 
1.36-10^^ 0. Besides the direct measurement 
of {ijjip) above, we may also determine {ipi/j) from 
the mf dependence on the pion mass, using Gell- 
Mann-Oakes-Renner (GMOR) relation: 



J TT 



48(m/ 



(1) 



The values obtained using both methods are given 
in lattice units in Tab. |21 Apart from the 3GeV 
case, the values from each method agree nicely. 

3. Non-perturbative renormalization 

In order to get physical values for continuum 
observables, results obtained from lattice simu- 
lation need to be renormalized. The RBC group 
has used the non-perturbative RI renormalization 
scheme which benefits from the 0{a) off-shell im- 
provement of DWF. A detailed explanation of this 
method is given in the paper 0]. There are two 
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Figure 1. Zs/Zg and Zp/Zg plotted versus for a ^ = 3GcV(lcft) and dynamical fermions (right). 
Table 2 

Values of (i/'V') from the GMOR relation, jof^ and from direct measurement (the last column). 



Action 



ruf 



DBW2 quench DWF 2GeV m/ > 0.01 

DBW2 quench DWF l.SGeV to/ > 0.02 

DBW2 quench DWF 3GeV to/ > 0.008 

DBW2 dynamical DWF 1.75GeV tod > 0.02 

systematic errors in this approach: lattice arti- 
facts which appear at large momenta and non- 
perturbative phenomena which appear at small 
momenta. I will describe each of these in turn. 

This approach gives renormalization fac- 
tors for quark bilinear operators uV id, {Ti = 
{1, 7i,, 75, 7,^75, cr^i/}) where u and d are quark 
fields, using a corresponding vertex amplitude A^. 
The renormalization factor Zi is determined by 
the condition: 



ZviZg ^Ar,,o(p,p)|p2 



1. 



(2) 



For the case of A5 and Ap, there are non- 
perturbative contributions which must be re- 
moved. They are shown in the equations: 



h.p,iau{ap,ap) 



^s.iatt{ap,ap) 



a^{qq) 



{apY{mf + TOre 



-ClZg 



+ ZmZq 



(3) 



CiZgda^qq) 
(ap)'^ dnrif 



Z,nZg. (4) 



Subtracting the effects of the first term in each 
of these equations, gives values of Zs/Zp close to 
unity over a large range of momentum, see Fig.^ 
Generally, renormalization factors in pertur- 
bation theory have logarithmic momentum de- 
pendence. This physical momentum dependence 
should be removed before attempting to identify 



2.3(2) nir^^ 2.0(3) • 10"* 

1.00(2) • 10-3 9.87(5) • 10"'' 

9.85(13) • 10-5 1.78(2) lO^* 

5.12(91) ■ 10-"^ 5.56(17) • 10-" 

errors. Therefore we define a scale invari- 
ant (SI) version of the RI renormalization factor: 

k^rl{{apf)^KrA{ap?)/CTA{ap?) (5) 



2.58(2) 
5.02(41) 
1.83(3) 
3.70(10) 

0(aV 



where Cr^ is calculated through three loops and 
is normalized so that Cviip^) = 1- As an ex- 
ample. Fig. 121 shows the bare Zs/Zg and its SI 
version for 3GeV and dynamical simulations. In 
this figure, the dashed line is a linear fit of Ap:'^ as 
a function of (ap)^, permitting us to remove this 
large momentum lattice artifact. 

To determine the conventional physical quanti- 
ties defined in the continuum theory, it is neces- 
sary to convert renormalization factors obtained 
in Rl-scheme to those in the MS scheme using a 
relation of the form: 



7RI 



1 



47r ° 



O^s y{2)RI 
\2^0 



ZMS ■ 47r " ■ (471)2"*^ ' ■ 

In order to extract the needed factor 
Zs and Zp from Eq. [2 we must determine Zq. 
This is best done indirectly by using ZA/Zg from 
the vertex function of the local axial current and 
Za obtained from comparing the local and the 
conserved axial current. This method is described 
in detail in Ref. ^ and gives values for Za of 
0.8876(3) for a'^ = 3GeV and 0.7576(6) the dy- 
namical case P|. Using these Za value, ZA/Zg 
and Zs /Zq from Rl-scheme and the Eq. |21 we de- 
termine the Z factors given in Tab. 0] 
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Figure 2. Bare and SI values Zs/Zq plotted versus for a ^ = 3GeV(left) and dynamical fermions (right). 
Table 3 

Results for (■0V') expressed in physical units 



Action 




m 


res 


Z{MS) 






DBW2 quench 


1.31(4) 


5.7- 


10-4 


0.699(16) 


0.265(6) GeV 


0.266(7)GeV 


DBW2 quench 


1.97(4) 


1.7- 


10-^ 


0.712(14) 


0.235(13)GeV 


0.246(9)GeV 


DBW2 quench 


2.86(9) 


9.72 


10-5 


0.832(8) 


0.345(ll)GeV 


0.282(1) GeV 


DBW2 dynamical mo = 0.02 


1.80(7) 


1.358 


• 10-3 


0.469(6) 


0.263(10) GeV 


0.253(18)GeV 



4. (V'V') values in physical units 

Using the renormalization factors discussed in 
the previous section, we obtain the values for 
(■)/)■!/') in physical units given in Tab. |31 The re- 
sults for (V'V') agree very well except for the direct 
calculation in the 3GeV case which is shifted 30% 
above the others. It is expected because of the a 
dependence of the explicit chiral symmetry break- 
ing in the domain wall scheme where can be 
written in physical unit as, 



coA^ 



Cl 



C2 



QCD 



2 „ — Lsct 



(7) 



The last term in this equation can be estimated 
to be of the same size as ^m^es- For the GeV 
case, this contributes an error to the value of 
(12(-0-0))^/3 quoted in Tab.Olon the order of 20%. 



Table 4 

Quenched 3GeV and dynamical Z factors, 
dynamical results are given for v/f' = 0.02. 



The 



Quantity 



3GeV quench 1.75GeV dyn. 



RI&SI ZA/Zq 

RIfcSI Zs/Zq 

Ms ZA/Zg 

MS Zs/Zq 



0.962(9) 
0.808(3) 
0.964(9) 
0.904(3) 



0.916(6) 
0.615(19) 

0.921(6) 
0.7366(98) 



This error is expected to decrease as — * cjo. 

5. Conclusion 

Values for {ipip) are obtained from different 
methods with different lattice scales. Good agree- 
ment is found in all but one case where a larger 
Ls is required. 

6. Acknov^rledgment 

I thank the RBC members for providing 
the simulation results analyzed in this re- 
port and many useful suggestions, especially 
Dr. C.Dawson for his help for NPR calculation, 
and Prof.N. Christ for his help. 

REFERENCES 



1. T. Blum, et al. RBC, hep-lat/0007038 

2. Contributions of Y. Aoki, C. Dawson, J. 
Noaki and K. Orginos, to Lattice Conference 
'02. Nucl.Phys.B.(Proc.Suppl) Vol.119. 

3. Contributions of Y. Aoki, C. Dawson, J. 
Noaki and K. Orginos, to Lattice Conference 
'03. 

4. T. Blum, et al. i^gg, .hep-lat/0102005| 



